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Abstract 

We consider the phase transition of a string with tension £\ to a string with a smaller 
tension £2 at finite temperature. For sufficiently small temperatures the transition 
proceeds through thermally catalyzed quantum tunneling, and we calculate in arbitrary 
number of dimensions the thermal catalysis factor. At £2 = the found formula for 
the decay rate also describes a break up of a metastable string into two pieces. 



1 Introduction 



A string-like configuration with string tension S\ can be metastable with respect to either 
a complete breaking [1, 2] or a transition to a string of a lower tension e 2 , with e 2 < £\ [3]. 
The metastability arises through existence of an energy barrier due to the emergence in 
the process of two interfaces (the ends of the original string) each having mass /i. At zero 
temperature the breaking of the string proceeds due to quantum tunneling and nucleation 
of a critical gap in the original string with the length £ c = 2fi/(ei — e 2 ), at which length 
the energy gain due to the lower tension in the gap compensates for the energy 2/x required 
for the formation of the interfaces. Once nucleated, the critical gap expands and converts 
the whole length of the original string into either a string with the lower tension, or into 
'nothing', which formally corresponds to the limit e 2 — > 0. This mechanism of the string 
transition bears strong similarity [1, 2, 3] to the decay of false vacuum[4, 5] and to the 
Schwinger process [6] of creation of pairs of charged particles by electric field. However there 
is an important difference between the string transition and the latter processes. Namely, 
the transverse waves on the string correspond to a presence of massless degrees of freedom, 
which substantially affect the probability of nucleation of the critical gap. The effect of the 
transverse degrees of freedom has recently been calculated [10] in terms of the preexponential 
factor in the rate of the nucleation at zero temperature of critical gaps per length of the 
original string in a rf-dimensional theory: 
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£i - £2. 



with fi R being the renormalized value of the mass of the interface that includes the effects 
of the transverse motion of the adjacent part of the string, and F is the factor contributed 
in the rate by each of the (d — 2) transverse dimensions and given by 




In this paper we consider the thermal effect in the rate of the string transition at low 
temperatures. This effect further exposes the difference between the breaking of a string 
on one side and the false vacuum decay and the Schwinger process on the other. Namely, 
as long as the temperature T is lower than the inverse of the critical length, T < l/£ c the 
thermal effects in the latter processes is exponentially suppressed as exp(— m/T) with m 
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being the lowest scale for particle masses in the theory, and these effects are essentially due 
to the (exponentially small) presence of massive particles in the thermal equilibrium[ll]. 
In the case of a string, however, the transverse waves on the string are massless so that 
their excitation has no suppression by the mass at arbitrarily low temperature. The thermal 
excitations of these waves create fluctuations in the distribution of energy in the string which 
catalyze the nucleation of the critical gap. Clearly, at T l/£ c the typical wavelength of the 
thermal waves A ~ 1/T is large in comparison with the critical length £ c , and the thermal 
effect in the rate is quite small, although not exponentially small. We find, as a result of the 
calculation to be presented in this paper, that the leading low temperature correction in the 
nucleation rate is given by the thermal catalysis factor 

We also find that as T approaches l/£ c the catalysis factor develops a singularity at £ C T = 1. 
At still higher temperatures the considered here string transition behaves, in a sense, simi- 
larly to the false vacuum decay [11], namely the regime of the transition changes to a different 
tunneling trajectory, so that the temperature dependence appears in the semiclassical expo- 
nential factor in Eq.(l), rather than in the preexponential term. 

The material in this paper is organized as follows. In the next Section we formulate 
the problem in terms of periodic configurations in Euclidean space-time, and in Sect. 3 we 
calculate the relevant path integrals and find the general formula for the catalysis factor in 
terms of expansion in powers of £ C T. In Sect. 4 we further analyze the general formula and 
find the first terms in the low temperature expansion of the catalysis factor, and also calculate 
this factor numerically for the temperatures approaching the critical point at £ C T = 1. 
Finally, Sect. 5 contains a summary of our results. 

2 Euclidean-space calculation 

The tunneling trajectory describing the decay of a metastable state can be found as a classical 
configuration, called "the bounce", in the Euclidean space-time [5], and the exponential factor 
for the decay rate is given by the classical action on the bounce. The fields in the bounce 
configuration approach their values in the metastable state at the boundaries of the space- 
time and they approach the final state values inside the bounce. 
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In the discussed here problem of the string transition the bounce is found in terms of the 
effective low-energy Nambu-Goto action for the states of the string. For two strings with 
tensions S\ and e 2 and a particle with mass \i propagating along the interface of the strings 
one can write the action in the form 

S = fiP + e 1 A 1 + e 2 A 2 , (4) 

with P being the perimeter of the interface and being the corresponding string's world 
surface area. The approximation of the action of the system with the effective low-energy 
one is applicable as long as massive internal degrees of freedom of the strings are not excited. 
This implies in particular that the approximation is applicable only as long as the thickness of 
the strings can be neglected comparing to other length scales in the problem, or equivalently 
only as long as all the momenta in the problem are small compared to the inverse thickness 
of the string. If the thickness of the string is r and the mass scale associated with it is 
M = l/r then we can write the condition of applicability of the approach presented in the 
form 

I > r , k < M (5) 

We consider initially a string with tension e 1 with very large length X stretched along the 
x direction. At zero temperature the stationary configurations for the action (4) are the 
trivial one, which is a flat world surface of the string lying in (t, x) plane, and the bounce 
configuration, which has a disk filled with the lower phase of the string with the radius 

R=-^, (6) 
£i - ^2 

as shown in Fig. 1. One can readily notice that the length of the critical gap is the diameter 
of the disk, £ c = 2R. The condition (5) for applicability of the effective action (4) requires 
that 

— = M R= ^— > 1 . 7) 

The rate of the critical gap nucleation is then found[7, 5, 8] by calculating the path 
integral Z\ 2 around the bounce and comparing it with the path integral Z\ around the 
trivial configuration: 

where A is the total area of the world surface for the string. It can also be reminded that, as 
explained in great detail in Ref . [8] , the imaginary part of Z 12 arises from one negative mode 
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Figure 1: Bounce configuration for zero temperature 



at the bounce configuration, and that due to two translational zero modes the numerator 
in Eq.(8) is proportional to the total space time area A in the (t,x) plane occupied by the 
string, so that the finite quantity is the transition probability per unit time (the rate) and 
per unit length of the string. 

The formula in Eq.(8) corresponds to a calculation of the decay rate as the imaginary 
part of the energy of the initial string. At a finite temperature T the corresponding relevant 
quantity is the imaginary part of the free energy [9], which one can calculate in the Euclidean 
space by considering periodic in time configurations with the period (3 = T~ l . In other 
words the thermal calculation corresponds to the path integration in the Euclidean space- 
time having the topology of a cylinder. The nucleation rate is then described by the same 
formula (8) with the action and the area being calculated over one period, A = X (3. 

In the present paper we consider only sufficiently small temperatures (3 > l c = 2 R, at 
which temperatures we find the thermal effects behaving as powers of T, which distinguish 
the string process from the decay of metastable vacuum[ll]. We also treat the length X 
of the metastable string as the largest length parameter in the problem, so that /J C I. 
Under these conditions the bounce corresponding to the action (4) is the same as at zero 
temperature, except that it is placed on a cylinder rather than on a large flat plane (Fig. 2). 

We aim at calculating the path integral over the variations of the string around the 
bounce configuration, which involves in particular the integration over the shifts of the 
string in transverse directions, i.e. those perpendicular to x. The contribution from each 
of the d — 2 transverse dimensions factorizes, so that it is sufficient to consider only one 
transverse dimension with further straightforward generalization to an arbitrary dimension. 
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Figure 2: Bounce configuration for nonzero temperature 



The coordinates on the cylinder (or on the plane where all the points separated by n/3 along 
Euclidean time are identified) are t, x, with t being the periodic time coordinate, and the 
coordinate orthogonal to the surface of the cylinder is z. The boundary conditions for the 
configurations over which we integrate are 

z(t,x = ±—^J=0, z(t + (3, x) = z(t, x) . (9) 

Introducing the polar coordinates (r, 6) in the (t, a:)-plane one can consider small variations 
of the classical configuration (Figs. 2, 3) in the following form. Variations of the interface 
are given by 

r(6)=R+f(6) z{R,9) = t{e), (10) 

while the variations of the surfaces of the string in the bulk are Zi(r, 8) and z-i (r, 9) with the 
boundary conditions on the interface 

z 1 (R,9) = z 2 (R,6)=m- (11) 

One can certainly notice that the symmetry of the polar coordinates is that of the disk 
describing the bounce, but it is not the symmetry of the periodicity condition in t. It will 
be clear from the subsequent calculation, that in fact all the discussed thermal effect can be 
viewed as originating from this mismatch, once it is properly accounted for. 
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Figure 3: Variation of the bounce configuration 



fn terms of the introduced variables the action (4) can be written in the quadratic ap- 
proximation in the deviations from the bounce as 

e x - s 2 2 J V ' 

where the primed and dotted symbols stand for the derivatives with respect to r and 9 
correspondingly, and the "out" and "in" integrals are taken over the surface respectively 
outside and inside the disk r < R. 

Similarly, the action around a flat initial string configuration in the quadratic approxi- 
mation takes the form 

S 1 = £ 1 Xf3 + ^ J rdrdO (V 2 + , (13) 

with the integral running over the whole space-time cylinder and z(r, 6) parametrizing small 
deviations of the string in the transverse direction. 
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3 Integration 



The path integral Z\ runs all transverse fluctuations z(r, 6) with the action given by Eq.(13). 
This integral contains no interface at r = R. However, as explained in Ref. [10], for the 
purpose of calculation of the ratio of the partition functions in Eq.(8) it is helpful to organize 
the integration for Z\ in the following way. First one fixes the values of z on the circle of radius 
R: z(R,9) = and integrates over the remaining bulk variables with fixed £(0). The 

full partition function Z\ is then found after path integration over the boundary variations 
£(#). As demonstrated in Ref. [10] this procedure in fact factorizes the full path integral 
into a product of "bulk" and "boundary" factors with the bulk factors being the same in 
Z l2 and Z\. As a result the ratio of the full path integrals in Eq.(8) is determined only by 
the boundary terms: 

dT = £i - gg ^ / 7T/i 2 \ ^(boundary) ^ 
di 27T y E 1 — E 2 J ^l(boundary) 

with the boundary partition functions 2i2(boundar y ) and ^(boundary) being given by 

-2l2(boundary) = j V£ 6Xp - — j dB ? - R j dQ {^d r Z 2c \ r=R - ^d r Z lc \ r=R ^j £ 

^i(boundary) = j T>£ exp -R y J dQ (d r z 2c \ r=R - d r z lc \ r=R ) £ , (15) 

where the functions z lc and z 2c satisfy the Laplace equation Az ic = with the boundary 
conditions 

z lc (R,9) = z 2c (R,9)=m . (16) 

The outer solution z\ c is also periodic in time z\ c {t + f3,x) = zi c (t,x) and satisfies the zero 
boundary condition at the spatial infinity 

z lc (t,x = ±—)=0 (17) 



2 J 

while the inner solution z 2c is required to be regular inside the disk. 

In order to do the path integrals one can expand the boundary function £(#) in angular 
harmonics: 

a i oo 

£(0) = -£= + -= £ [a n cos(n6) + b n sin(n0)] . (18) 

For the inner solution z 2c to the Laplace equation, i.e. at r < R one finds no difficulty in 
finding the harmonics matching the boundary function at the interface [10]: 

4? = , 4? (r, 6) = 4= K cos(n0) + b n sin(n0)] £ . (19) 
V27T V 71 " 
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For the outer solution z\ c however there is a difficulty due to the previously mentioned 
mismatch between the symmetry of the boundary and of the periodicity conditions. It is 
impossible to choose the solution to the Laplace equation for outer string bulk variable to be 
z ic( r > 6) ~ r_r \ since it is not periodic in time. In this situation in order to have a periodic 
solution one can perform a periodic mapping of the cylinder on the plane and consider the 
outer solution of the Laplace equation as the sum of the solutions produced by a "source" 
at each period as illustrated in Fig. 4. Introducing the complex variable w = t + % x, we 
construct the solutions for the functions z\ c (r, 9) using the harmonic real and imaginary 
parts of the following basis set of periodic functions, satisfying the boundary condition (17) 
at large \x\, 




Figure 4: Periodic configuration on the plane 

Clearly these functions are periodic in the Euclidean time by construction, with the period 
(3. Also the functions gk{w) with k > 1 are analytic complex functions, so that their real 
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and imaginary parts are harmonic. In the functions go{w) and gi(w) the explicit dependence 
on x, introducing non-analyticity, is linear and is thus also harmonic, so that their real and 
imaginary parts do satisfy the Laplace equation. 

An arbitrary periodic outer solution to the Laplace equation, satisfying the boundary 
condition (17) at large \x\ can be expanded in the series 

oo 

z lc (r, 6) = A Reg (w) + A k Reg k (w) + B k lmg k (w), (21) 

k=i 

The disadvantage of this set of solutions (hence of such an expansion) is that it is not 
orthogonal, so that the expression for the action up to quadratic terms is not diagonal in 
this basis. Therefore, the calculation of the integral is not just a calculation of the product 
of eigenvalues. To find the integral over the amplitudes of the Fourier harmonics one has to 
express the amplitudes of the solutions chosen A k ,B k in terms of the amplitudes a k , b k . The 
relation between the coefficients A k ,B k and a k , b k can be found from the matching condition 
(16) on the interface 

— + — J2 K cos(n6>) + b n sin(n0)] = A g (w) D + [A k Reg k (w) + B k lmg k (w)] 
v2tt V 7r „=i r=R fc= i r=R 

(22) 

One can readily notice that the function g contains a large constant term, proportional to 
X, which totally dominates the matching condition for the ao mode at r = R, so that 

-A ^[l + O(R/X)}. (23) 



v 7 ^ u 2(3 

For this reason the effect of the mixing between ao and higher modes is suppressed by inverse 
powers of X and can be ignored in the limit of a long string. For this reason in considering 
the mixing of the modes in the following calculation we keep only k ^ 0. Furthermore the 
linear in x terms in the functions g and g\ are suppressed at r = R by the factor R/X and 
we also neglect them. 

In what follows we consider the expansion of the functions g k at r = R in powers of Rj (3, 
which expansion, as will be seen later, converges at R < (3/2. Using 

oo 

(l + ^-^l + Et-ir^C^, (24) 

p=l 

where Cp +fc _ 1 are binomial coefficients, 

p! (jfc_i)! W 
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and also a definition of the Riemann (^-function 



c(*) = E»" 



we find for k ^ 



with 



d 



pk 



'a 



(26) 



(27) 



(28) 



We have omitted in the expression (27) a constant term, which describes the mixing with 
the do mode as well as a term explicitly proportional to R/X. Using the expansion (27) for 
gk{w) and considering the real part of the functions g k (w), one can find the coefficients a; in 
terms of A k 

ai = - V A k / Reg k cos(W)dO = A x + V A k d ik , 



or in the matrix form 

A=(l + D)~ 1 a , 

where matrix D has elements d\ k . Similarly for the imaginary part one gets 

= r lmg k cos(W)d9 = Bt-J2 B k d L 

and 

B = (1 - 6 . 



(29) 
(30) 



Hk 



(31) 



(32) 



As usual, the contributions from cos and sin modes are independent, and we consider the 
contribution to the boundary term in (15) from the even (cos) harmonics first 
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(Rg^d r 



R 



d6 



t-2-K 



d6 



E a « cos(W) 



E 

k 



k[ — J cos(k9) -J^pdpk 
P =i 



k 



ka k - E pa P d p k 



Introducing the matrix 



k, P 



A> = diag(l,2,...,n,...) , 



(33) 



(34) 
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one can rewrite the expression (33) in the matrix form 



r2n 

-J (Rg {R) d r g iR) )\ R d6 = aN{l- D)A . 
A substitution in this expression of the solution for A in terms of a (30) leads to 
-J (Rg (R) d r g (R) )\ R d9 = aN(l - D)(l + D)~ l a . 



(35) 



(36) 



Clearly, for the odd (sin) modes one gets the same expression with the replacement D — > —D. 
Collecting all the terms together one can write the result for the boundary partition functions 
(15) as 



-^(boundary) — Det 



Z\ (boundary) — Det 



(ei - e 2 )N 2 + e 2 N + £l N(l - D) 



1 



l + D 



-1/2 



{D —D} 



£lJV + £iJV(l - D) 



-1/2 



{D^-D}, 



(37) 



1 + D\ 

where {D — > — D} means that one should take the expression and make the replacement 
D -> —D. 

The zero temperature limit for the probability rate dT /d£ formally corresponds to setting 
D — > 0. Thus one can use the known result for the zero temperature decay rate [10], and 
concentrate on a calculation of the thermal catalysis factor K defined as 



dl dl 



(38) 



In a d- dimensional theory, i.e. with d — 2 transverse dimensions, the catalysis factor can be 
written as K = G d ~ 2 , where G is the factor per each transverse direction given by 



G = 



7 Z D=0 
^-'12(boundary) l(boundary) 

;rD=0 "t 7 " " 

^12(boundary) -^l(boundary) 



(39) 



According to Eq.(37) it is a matter of simple algebra to express the factor G in terms of the 
matrix D: 

( N — 1 



G = Det 

where we have introduced the parameter b 



N + b 



D 



-1/2 



(40) 
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4 Analysis of the general formula 



In this section we consider in some detail the temperature effect in the string transition rate 
described by our general formula (40) in the situation where the inverse temperature is larger 
then the diameter of the classical configuration (bounce) (3 > 2R. We first notice that due to 
the presence of the factor (JV — 1) the first elements from the first row and the first column 
of the matrix D, d lk and d p i, enter the expression [(N — 1) D] 2 with zero coefficients, so that 
the final result (40) does not depend on them. Furthermore, one can see from Eq.(28) that 
the matrix element d p k is not equal to zero only if the indices p and k have the same parity. 
Hence, there is no mixing between the amplitudes with even (a 2 i, b 2 i) and odd (a 2 i+i, &22+1) 
indices. Therefore the determinant in the (40) can be written as a product of determinants 
corresponding to even and odd amplitudes 



Dct 



1 - 



N-l 
N + b 



D 



Dct 



1 - 



'2N-1 
2N + b 



U 



Det 



1 - 



2N 



21 



2N+l+b 



V 



(41) 



where the matrix elements of the matrices U and V are 



Upk — d 



2p2k , 



V pk = d 



2p+l,2fe+l , 



(42) 



and the indices p and k take values 1,2,.... 

For practical calculations it is also convenient to write the expressions for the elements 
of the matrices entering in Eq.(41) in terms of their indices: 

2p+2k 



2N-1 
\2N + b 
2N 
2N+l+b 



U 



V 



pk 



pk 



2p-l (2 P + 2fc-l)! /_ 
2p + b(2p)\(2k-l)\ a2p + 2k} {(3) 

Ap (2p + 2k + l)\ (R 
2 P +l + b(2p+l)\(2k)\ a2p + 2k + 2} {(3, 



2p+2k+2 



,(43) 



with p, k — 1, 2, 

In order to find the first thermal correction at low temperature one can expand the 
matrices in power series using well known formula for the determinant 



Det(l — M) — exp [Tr In (1 - M)\ = exp 



00 M i 



1=1 



l-TrM + 0(M 2 ). (44) 



12 



In our case 



2N- 1 
2N + b 



U 



I 36C 2 (4) (R\ 600 C(6) 2 ( R\ 120C(4)C(6) (R 
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(2 + 6) 2 \pl (2 + 6) (4 + 6) \ p 



180C(4)C(6) (R 
(2 + b) (4 + b) U 
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(2 + 6) 2 \PJ 
600 C(6) 2 (R 



12 



(2 + 6) (4 + 6) 1/3 



2iV 



2iV + 1 + 6 



2 / 1600 C(6) 2 /'iZ' 
I = 

V 



12 



(3 + 6)2 \p 



■ J 

(45) 

(46) 
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Therefore the first correction to the zero temperature value of the rate is proportional to 
R 8 /P 8 and is given by Eq.(3). 



1.4 
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1.2 
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-ft— 8- 



0.35 0.40 0.45 



0.50 
R T 



Figure 5: The thermal catalysis factor per each transverse dimension vs RT for 6 = 1 (circles) 
and for 6 = 10 (solid circles) 



The series for the function G(T) diverges when (3 = 2 R. The corrections for the tem- 
perature close to (2R)~ 1 can be found numerically. The proximity to this point defines the 
number of terms which should be taken into account in the series for G(T). The plot for the 
function G(T) vs the parameter R//3 = RT calculated numerically with 50 first rows and 
columns retained in the matrices U and V is shown in Fig. 5. 
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5 Summary 

In the present paper we considered the transition between two states of a string at finite 
temperature. The result for the rate of the transition in arbitrary number of dimensions is 
obtained in a form of a product of the value at zero temperature and the catalysis factor (40). 
The numerical analysis of the expression shows (see Fig.4) that for smaller b the influence 
of the temperature on the process becomes noticeable earlier and the catalysis factor grows 
more rapidly with temperature. Our result is applicable in the range of temperatures for 
which the associated period in Euclidean time is larger then the critical size of the classical 
configuration £ C T > 1, so that it is still possible to inscribe the bounce in one period. We 
have shown that, similarly to the process of false vacuum decay at finite temperature in two 
dimensions [11], there is no change in the exponential behavior of the probability rate while 
the temperature is small enough, and that such change starts only when (5 becomes smaller 
than the diameter of the bounce. On the other hand we have demonstrated that unlike 
in the case of false vacuum decay, there arise thermal corrections to the transition rate at 
low temperature, which corrections have power dependence on the temperature and are due 
to thermal waves excited on the string at arbitrarily low temperatures. In practical terms, 
as can be seen from Fig. 5, the value of the catalysis factor changes from 1 for RT = to 
approximately 1.5 for RT = 0.49, i.e. quite slowly, even for b = 1. However there is an actual 
singularity in this factor at RT = 1/2 so that the string transition becomes significantly 
accelerated by thermal effects for temperatures close to l/£ c . 
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